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Neutralino Annihilation into Massive Quarks with SUSY-QCD Corrections 
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We compute the full 0(a„) supersymmetric (SUSY) QCD corrections for neutralino annihilation 
into massive quarks through gauge or Higgs bosons and squarks in the Minimal Supersymmetric 
Standard Model (MSSM), including the known resummation of logarithmically enhanced terms. The 
numerical impact of the corrections on the extraction of SUSY mass parameters from cosmological 
data is analyzed for gravity-mediated SUSY breaking scenarios and shown to be sizable, so that 
these corrections must be included in common analysis tools. 



INTRODUCTION 

The astrophysical and cosmological evidence gathered 
in recent years supports the hypothesis of dark mat- 
ter in our Universe. In particular, models of structure 
formation favor Cold Dark Matter (CDM) consisting of 
Weakly Interacting Massive Particles (WIMPs) with non- 
relativistic velocities. The five-year data of the Wilkin- 
son Microwave Anisotropy Probe (WMAP), combined 
with the results of supernovae experiments and baryonic 
acoustic oscillation data, constrain the relic density fl of 
CDM in the Universe at 95% (2a) confidence level to [1| 



0.1097 < n 



CDM 



h 2 < 0.1165. 



(1) 



Here, h denotes the present Hubble expansion rate Hq in 
units of 100 km s _1 Mpc -1 . 

In contrast to the Standard Model (SM) itself, its 
various extensions can provide viable candidates for 
WIMPs. In the Minimal Supersymmetric Standard 
Model (MSSM) with i?-parity conservation, the Light- 
est Supersymmetric Particle (LSP) is such a candidate, 
if it is a color singlet and electrically neutral. In many 
scenarios, in particular those where SUSY breaking is 
mediated by gravity, the LSP is the lightest neutralino 
and thus a suitable dark matter candidate in a large part 
of the parameter space. To calculate the number density 
n of the relic particle, one has to solve the Boltzmann 
equation 
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where n oq is the density of the relic particle in thermal 
equilibrium and v is the relative velocity of the annihi- 
lating pair. 

The thermally averaged annihilation cross section 
(&a.nnv) includes all (co-) annihilation processes of the 
dark matter particle into SM particles. It is dominated 
by two-particle final states, most notably by fermion- 
antifermion pairs and by combinations of gauge (W^, Z°) 



and Higgs bosons (h°, H°, A°, if*) 2\. The nature of the 
final state contributing most depends strongly on the re- 
gion of parameter space. Fermion final states have the 
clear advantage that they are always kinematically al- 
lowed. Their leading contribution is proportional to the 
mass of the fermion, so that we focus our attention here 
on the annihilation into the massive quarks of the third 
generation. 

The relic density of the lightest neutralino depends on 
the SUSY-breaking parameters of the MSSM, which de- 
termine the nature of the neutralino as well as the cou- 
plings and masses that appear in the annihilation cross 
section. The number of free parameters is often reduced 
to a few universal parameters imposed at the high-scale 
e.g. the five parameters m , m l / 2l A$, tan/?, and sgn(/i) 
in minimal supergravity (mSUGRA). Using the experi- 
mental limits in Eq. |T]), one can then constrain the pa- 
rameter space in a complementary way to collider and 
low-energy experiments. 

This analysis is made possible by public computer 
codes, which perform the calculation of the dark mat- 
ter relic density within models beyond the SM. The 
most popular and developed codes are DarkSUSY 0] and 
micrOMEGAs |4]. The implemented processes are mostly 
calculated at leading order, and higher order corrections 
are only included for some very sensitive quantities. How- 
ever, owing to the large magnitude of the strong coupling 
constant, all QCD and SUSY-QCD corrections signifi- 
cantly affect the annihilation cross section. 

In this Letter, we extend our previous work 

i by com- 
puting the full QCD and SUSY-QCD corrections to neu- 
tralino annihilation into top and bottom quark- antiquark 
pairs. As an example, we present the impact of the radia- 
tive corrections in mSUGRA scenarios where the annihi- 
lation into top quark pairs is dominant. Moreover, we in- 
clude a study of bottom quark final states at small tan /?, 
where we find agreement with previous results of a sim- 
ilar calculation 6|. Since light quark final states do not 
lead to sizable contributions in the analyzed mSUGRA 
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FIG. 1: Tree- level diagrams for the annihilation of neutralinos 
into quarks. The crossed u-channel diagram is not shown. We 
denote the neutral Higgs bosons (h°, H°, A ) by H®. 
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5800 


536 


-1500 


10 
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0.111 




72% 
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3200 


520 





50 


+ 


0.110 


8% 


71% 



TABLE I: High scale mSUGRA parameters (in GeV) together 
with the corresponding neutralino relic density and the con- 
tributions of the massive quark anti-quark final states to the 
annihilation cross section for our selected scenarios. 




FIG. 2: SUSY-QCD loop diagrams contributing to the anni- 
hilation of neutralinos into quarks. We leave out the diagrams 
for the M-channel contributions and the crossed box diagrams. 



scenarios, their discussion is postponed to a later publi- 
cation. 



CALCULATION DETAILS 

The annihilation of neutralinos into quarks proceeds at 
tree-level through an exchange of the Z-boson and Higgs- 
bosons in the s-channel or of scalar quarks (squarks) in 
the t- and it-channels (see Fig. [1]). To reach a sufficient 
annihilation rate, the cross section has to be enhanced 
either by an s-channel resonance or by a small squark 
mass. 

For these processes, we calculate the full one-loop QCD 
and SUSY-QCD corrections involving the exchange of a 
gluon or a gluino. These affect the vertices including Z- 
and Higgs-bosons and a quark-antiquark pair as well as 
those with neutralinos, quarks, and squarks. In addition, 
one has to include the box diagrams arising from the ex- 
change of a gluon or a gluino between the final state 
quarks in the t- and u-channels (see Fig. [2]). All relevant 
divergent integrals are evaluated in the modified dimen- 
sional reduction scheme (DR) . We use on-shell renormal- 
ization of the wave functions and couplings in order to 
eliminate the ultraviolet divergencies. The squark mix- 



ing matrix is also renormalized on-shell as proposed in 
Ref. Q. In the Yukawa couplings, however, we use the 
quark masses defined in the DR renormalization scheme. 
Full details of the calculation will be presented in a forth- 
coming paper. 

For the bottom quark, we start from the input value 
(nib), evolve to the scale Q using three- loop SM 
renormalization group equations, change from MS to DR 
scheme using the corresponding relation at two loops [8(] , 
and finally include the MSSM threshold corrections com- 
prising the sbottom-gluino and stop-chargino one-loop 
contributions. The latter are considerably enhanced for 
large tan j3 or large At, so that they have to be resummed 
to all orders of perturbation theory . Denoting the re- 
summable part by A& and the finite one-loop remainder 
by Am&, the bottom quark mass is then given by 



MSSM 



(Q) = 



,DR 
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(3) 



The DR-mass of the top quark is obtained from the 
on-shell value m t — 172.4 GeV [k} by subtracting the 
finite part of the on-shell mass counterterm. 

The infrared divergence connected to the exchange 
of a massless gluon is also regularized dimensionally. 
To cancel the infrared divergent poles, we include the 
bremsstrahlung process with an additional gluon in the 
final state. To allow for integration of the divergent 
bremsstrahlung matrix element and for the cancellation 
of the infrared divergencies, we use the dipole subtraction 
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formalism for massive partons 

We have performed several checks of the calculation. In 
particular, we have compared our analytical calculation 
of the virtual part with the one produced by the packages 
FeynArts and FormCalc 12j and the QCD corrections 
to the s-channel Higgs-boson exchange with the known 
results of Ref. 0. 



NUMERICAL ANALYSIS 

To evaluate the impact of the QCD and SUSY-QCD 
corrections on the neutralino relic density, we have cho- 
sen three typical parameter points in mSUGRA shown in 
Tab. [U They were selected so that the top and bottom 
quark anti-quark final states dominate the total annihi- 
lation cross section and the relic density lies within the 
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FIG. 3: The effects of the radiative corrections on the cross 
section as a function of center-of-mass momentum (top), the 
prediction of the neutralino relic density as a function of the 
gaugino mass parameter miy 2 (center), and on the cosmo- 
logically favored regions in the mo-mx/2 plane (bottom) for 
our parameter point 1 with annihilation into bottom quarks 
through the light Higgs boson resonance. 



WMAP range of Eq. Q. Moreover, the SUSY particle 
mass spectrum resulting from the high scale parameters 
respects the experimental mass exclusion limits, and the 
branching ratio of the decay b — > sj lies within 2a of the 
experimental bound BR{b -> s 7 ) = (3.52 ± 0.25) • 1(T 4 , 
obtained from combined measurements by BaBar, Belle, 
and CLEO 



The high-scale parameters are evolved down to the 
electroweak scale using SPheno 2.2.3 [l5| . The neu- 
tralino relic density is calculated with micrOMEGAs 2. 1, 
where we have included the full radiative QCD and 
SUSY-QCD corrections to the annihilation cross section. 
Our first parameter point at low tan/3 has been chosen 
near the bulk region with a low fermion mass parameter 
nii/2, but a rather large scalar mass mo in order to avoid 
the constraint coming from BK(b — > s 7 ) . The remaining 
top-quark dominated points lie both in the focus point 
region, where mo is very large. At large tan/3, there is 
also a contribution from bottom quarks due to their im- 
portant coupling to the CP-odd Higgs boson. 

In the first panels of Figs. [SHU we show the numerical 
results for the relevant annihilation cross sections as a 
function of the relative momentum p cm , that is related 
to the center-of-mass energy -^/s and the neutralino mass 



through s = 4 m; 



+ m\) 



Shown are the leading or- 



der results with DR Yukawa couplings (dash-dotted), the 
approximation already included in micrOMEGAs (dashed), 
and the result including our full one-loop QCD and 
SUSY-QCD corrections (solid). We also show, in arbi- 
trary units, the Boltzmann distribution function involved 
in the calculation of the thermal average at the freeze-out 
temperature (shaded area). It indicates which centre-of- 
mass momenta contribute most to the relic density. E.g. 
in Fig. [3] one sees that an important contribution comes 
also from energies below the light Higgs boson resonance, 
which is clearly visible on the top panel. 

The center panels of Figs. [3H5] show the corresponding 
predictions for the neutralino relic density as a function of 
the gaugino mass parameter around our parameter 
points 1, 2, and 3, respectively. The horizontal bands 
indicate the experimental limits of Eq. ([1]). Due to the 
enhanced cross section, the relic density is reduced by 
about the same amount, and the region of parameter 
space agreeing with the constraints is modified. 

In the bottom panels we finally show the cosmolog- 
ically favored regions in the mo-m^ plane that re- 
sult from a relic density calculation using the cross sec- 
tion at the tree-level, the approximation included in 
micrOMEGAs, and the full one-loop SUSY-QCD corrected 
cross section. Note that the effect of the corrections is 
larger than the experimental errors, so that three dis- 
tinct bands are observed. The grey shaded regions corre- 
spond to the points that are excluded due to mass limits, 
no electroweak symmetry breaking (EWSB), or the con- 
straint from BR(6 — > sj). Concerning the "bulk" region 
at low tan/3, the favored band is shifted to smaller val- 
ues of rrii/2 and thus towards the region excluded by LEP 
mass limits. The cosmologically favored focus point re- 
gion, however, is found at higher values of m.1/2 and lower 
values of mo after including the SUSY-QCD corrections. 

Our additional corrections to the relic density for the 
case shown in Fig. [3] amount to about 10% with respect 
to the approximation implemented in micrOMEGAs. This 
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FIG. 4: Same as Fig. 1 for our scenario 2 with neutralino FIG. 5: Same as Fig. H for our scenario 3 with neutralino 
annihilation dominantly into top quarks. annihilation into top and bottom quarks. The corrections to 

the bb final state only are also shown separately. 



is due to the low value of tan (3, which suppresses the 
part of the usually dominant resummed correction which 
is included in micrOMEGAs, making the rest of the SUSY- 
QCD correction relevant. In the case of top quark final 
states, shown in Figs. H] and El the additional correc- 
tions can amount up to 15%. Note that SUSY-QCD 
corrections for top quark final states are not considered 
in micrOMEGAs. It is interesting that in the focus point 
scenario with large tan /3, the correction to the bottom 
quark final state has a sizable contribution to the total 
correction. As a result, the relative contribution of the 
annihilation into bottom quarks increases compared to 



the 8% at tree-level. This is well visible in the centre 
and bottom panels of Fig. [5] and shows the importance 
of the tan (3 enhanced corrections to the bottom Yukawa 
coupling. 

CONCLUSIONS 

In summary, we have calculated the full one-loop QCD 
and SUSY-QCD corrections to neutralino annihilation 
into third generation quarks. A numerical evaluation has 
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shown that the corrections have sizeable effects on the an- 
nihilation cross section and in consequence on the extrac- 
tion of SUSY mass parameters from cosmological data 
assuming that the lightest neutralino is the cold dark 
matter candidate. The induced difference is of the same 
order of magnitude as the experimental error from cosmo- 
logical precision measurements, so that the full one-loop 
corrections should be taken into account when analyzing 
the SUSY parameter space with respect to cosmological 
data. 
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